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In this paper, model-predictive control (MPC) is proposed for controlling power source of accelerators. The 
system state equation is employed as the predictive model. With MPC, the difference between possible output 
and the ideal output is forecasted and decreased, so that the system can trace the ideal trail as closely and 
quickly as possible. The results of simulations and experiments show that this method can reduce influence of 


low frequency noise. 
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I. INTRODUCTION 


Particle accelerators are widely used in physics research, 
clinic treatments and environmental protection. Performance 
of the power source is important to accelerator operation as it 
determines how well the magnetic field accelerates the parti- 
cles, hence the high performance requirement on DC power 
supply of accelerators [1, 2]. 


If the stability requirement is just 1073, a normal negative 
feedback control will be good enough. When a power stabil- 
ity of 10-°~10~° is required, however, the changes in am- 
plification, power grid noise, and many other factors, begin 
to affect the power stability. Against such a requirement on 
power stability, the PID (proportion, integration and differen- 
tiation) method is often used, because it is easy and effective. 
However, this method cannot satisfy requirements all at once 
on accuracy, quickness and controlling simplification [3-7]. 


Model-predictive control (MPC) is one of the modern con- 
trol theories developed based on the rapid progresses of com- 
puter technology. It grew quickly in 1970’s and has been 
widely used in industries. One of its advantages is that MPC 
does not need high precision model demand for the controlled 
system. By using the rolling optimization, the error of actual 
output from the ideal output is minimized. The system’s next 
prospective output will be justified in advance by the feedback 
of the error. Therefore, influence of the model uncertainty, 
undesirable noise and digital time delay can be reduced, and 
the system robustness is improved [8-13]. 


In this paper, together with the noise, the given input is 
added to the power source after it is modified with MPC. 
The control system is a close-loop structure. The results of 
a Matlab simulation show that influence of the noise at low 
frequency can be reduced effectively. A real model is set up 
and FPGA in VHDL language is devised so that the control 
method can actualize on it. Test results show that the method 
is effective. 
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Il. EQUIVALENT MODEL OF CURRENT POWER 
SOURCE 


Generally, a current power source can be simplified as in 
Fig. 1. L3, R3, i3 and Vo are the inductor, resistor, and 
output current and voltage, respectively, of the magnetic field 
coil; C1, Co, L1, Rı and Rə form the filter section. 


Fig. 1. Structure of the current power source. 


When T7} and T; is conducting, V; = E; and when To 
and T3 is conducting, V; = — E; so the following equivalent 
equations can be derived: 


Vi = Riit + Lidizi/dt, (1) 

in, = CdVo/dt + iz3 + C2dVo, /dt, (2) 
Vo = izR + Ladig/dt, (3) 

Vo = Vo, + RoC2dVo, /dt. (4) 


From Egs. (1)— (4), the state equation can be obtained: 
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where, V2 = E(t. — t1)/(t2 + tı), tı and tg are conducting 
time of T) and T4, and To and T3, respectively, within one 
time period. 


+ (0) V1, (6) 


I. ESTABLISHMENT OF PREDICTIVE CONTROL 
MODEL 


The overall MPC control structure is shown in Fig. 2. A 
predictive model is designed to forecast a prospective output 
Yp(n + 1). As the predictive model is not the real model, so 
the output yp(n + 1) is modified to get the final prospective 
output Ypr(n + 1), from the real model output y(n). Then an 
optimization is performed to minimize a cost function of the 
tracking error and the finale input u(n) is obtained. At every 
new sampling period, the measured real output modifies the 
forecasted output on the basis of the predictive model. Then, 
a new round of optimization follows, and this forms a close- 
loop control [13-20]. 


Feedback 
adjustment 


Fig. 2. Structure of MPC. 


A. The predictive control model 


The predictive control model forecasts the subsequent out- 
put of the system on basis of the actual output just measured 
and the future time inputted. Any kind of model having some 
features of the system can be used as predictive model. In 
this paper, the predictive model is state equation of the power 
source. The discrete state equations can be described as: 


a(n+1) = Ga(n) + Hu(n), (7) 


y(n +1) = Ca(n +1), (8) 


where G = I + AT, and H = BT. I is the identity matrix, 
T is the sample period. A, B, C is state matrix, input matrix 
and output matrix of the system’s state equation, respectively. 
Then, the predictive output in time n + 1 can be obtained: 


Yp(n + 1) = CGa(n) + CHu(n), (9) 


where x(n) and u(n) are state variables and input at n time. 
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B. Feedback adjustment and rolling optimization 


In the MPC, feedback adjustment and rolling optimization 
are required. Since the predictive model is not the real model, 
so its output Ypr(m + 1) need to be justified by the error e(n) 
of the actual output from the last predictive output 


Ypr(n + 1) = y(n + 1) + he(n), (10) 


e(n) = y(n) — y(n), (11) 


where y(n) and yp(n) are the actual output and predictive 
output of n time, respectively; and h is error correction coeffi- 
cient (h = 1). This feedback adjustment makes the predictive 
model more assimilate to the real model. 


C. The rolling optimization 


The purpose of predictive control is to make the system 
output trace the ideal output as quickly and closely as possi- 
ble. So the control objective is to minimize the sum of square 
errors of the prospective value from the given value. The in- 
put, however, should not be too large. Thus, the quadratic 
performance expression is: 


J = plypr(n +1) — y(n + 1)? +qu(n)?, A2 


where p and q are weighting coefficients. It shows that as the 
process goes on, the optimum varies on line. So this opti- 
mization relates only to dynamic performance of the system. 

To minimize the objective function, the derivative of J 
should be equal to zero and then the control value is: 


u(n) = zpCH/[p(CH)’ + q), (13) 
z=y(n+1) — CGz(n) — e(n). (14) 


IV. SIMULATION AND EXPERIMENTAL RESULTS 


Based on the model in Fig. 1, a simulation was performed 
with Matlab. Specifications of that model are: Lı = 0.3 mH, 
Rı = 0.01 Q, L3 = 91.4mH, Rz = 0.0796 Q, Cı = 10 uF, 
Cy = 47uF, and Rp = 1. The sampling frequency is 
10 kHz, and the total simulating time is 10s. As shown in 
Fig. 2, noise (n) of different frequencies was added to the 
power source to check their influence on the controlled ob- 
ject. The results are given in Figs. 3 and 4 

Six simulations were performed with two kinds of input: 
ramp input (Fig. 3) and sine input (Fig. 4). The noise added 
to the system was in sine wave. In order to find how dif- 
ferent noises affect the output, two kinds of frequency were 
included: high and low. The specific input and noise for each 
simulation are given in Table 1, where nj), stands for low 
frequency noise and npign, high. 
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Fig. 3. (Color online) Output of ramp input with (a) low frequency noise, (b) low frequency +30 Hz noises, and (c) low frequency +100 Hz 
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Fig. 4. (Color online) Output of sine input with low frequency noise. 


TABLE 1. The input and noise specification of the simulation 


Figs. Inputs Now Thigh 

3(a) u=t 1 x sin(27 * 2 x t) 0 

3(b) u=t 1x sin(2a*2x*t) 0.1 * sin(27 * 30 x t) 
3(c) u=t 1* sin(2m7 *2*t) 0.1 x sin(27 * 100 x t) 
4(a) u=10*sin(t) 3 *sin(27 » 2 x t) 0 

4(b) u= 10 xsin(t) 3*sin(2a7*2*t) 0.1 * sin(27 * 30 x t) 
4(c) u= 10 x sin(t) 3*sin(27*2%*t) 0.1 * sin(27 x 100 * t) 


In Figs. 3 and 4, the red curves are system output controlled 
by MPC, and the black curves are without the MPC. The X- 
axis is the time and the Y-axis is the system’s output of corre- 
sponding input. 

In Fig. 3(a), the input is in ramp wave. Compared to the 
red curve, the black one reflects the influence of the noise 
on the output. In fact the red curve is also the ideal output 
that has no disturbance on the system. So it is clear the MPC 
diminished the noise influence. In Fig. 3(b), in addition to 
the same input and noise as in Fig. 3(a), a 30 Hz noise was 
added to the system. We can see the red one is much more 
close to the ideal curve than the black one. This is because 


Fig. 5. The specification of the power source. 


that the influence of low frequency noise was greatly reduced 
by MPC. However, comparing to the ideal output in Fig. 3(a), 
there are ripple waves in the red curve, caused by the 30 Hz 
noise. In Fig. 3(c), the added is noise which is even higher 
(100 Hz), and bigger ripple waves are seen. So, the MPC 
decreases influence of low frequency noise greatly but it can 
do little against higher frequency noises. 

In Fig. 4, the input is sine wave. The red curve in Fig. 4(a) 
is the ideal output without the disturbance in the input. From 
Figs. 4(b) and Fig. 4(c) the red curves are much more similar 
to the ideal one than the black ones, but there are ripples, too, 
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Fig. 6. Experimental results without (a) and with (b) MPC. 


at noises of higher frequencies. The MPC works against low 
frequency noise of sine wave. 


A real model with specifications shown in Fig. 5 was made. 
Fig. 5 differs from Fig. 1 in that there is resistance for the 
capacity and inductance in the real model. Hardware imple- 
mentation that actualized the proposed method was achieved 
by using GX-SOC/SOPC-CIDE simulator box. The system 
input was sine waveform with some regular noise of differ- 


[1] Chen Y X, Yan H B, Huang Y Z. High Pow Laser Part Beams, 
2010, 9: 2138—2142. (in Chinese) 

[2] Xia J W, Zhan W L, Wei B W. Atom Energ Sci Technol, 2009, 
43: 150-158. (in Chinese) 

[3] Shi C F, Huang Y Z, Chen Y X. Atom Energ Sci Technol, 2013, 
47: 1652-1655. (in Chinese) 

[4] Wang B Y and Ma X L. Pow Electron, 2010, 2: 67—69. (in 
Chinese) 

[5] Wang R K, Huang Y Z, Chen Y X. Nucl Electron Detec Tech, 
2013, 33: 558—561. (in Chinese) 

[6] Huang Y Z, Chen Y X, Zhou Z Z. Nucl Phys Rev, 2011, 28: 
296-299. (in Chinese) 

[7] Wang B Y. Ph. D. Thesis, Tianjin University, 2006. (in Chinese) 

[8] Wang Z P, Xie Y X, Wang Y. Auto Inf Eng, 2013, 34: 34-40. 
(in Chinese) 

[9] Liu Y, Wang H Q, Li P. Cont Theory App, 2009, 26: 107-110. 
(in Chinese) 

[10] Morel F, Lin-Shi X, Rétif J M, et al. Electric Pow Syst Res, 
2008, 78: 1437-1447. 
[11] Grimble C J and Ordys A W. Annu Rev Control, 2001, 25: 


Nucl. Sci. Tech. 25, 050203 (2014) 


ent frequencies. Fig. 6 shows the experimental output wave- 
form without and with the MPC. One sees that the curve in 
Fig. 6(b) is more similar to the standard sine wave than that 
is in Fig. 6(a), because of the noise reduction by the MPC. 
Therefore, the proposed control strategy is validated. 


V. CONCLUSION 


A real object to be controlled is much more complicated 
than the theoretical model. Because of unexpected distur- 
bance, inaccuracy of components and computational delay, 
the real output differs with the ideal output. This problem 
must be taken into account in design of the controller algo- 
rithm to improve its performance and robustness. In this pa- 
per, we present an MPC method to reduce errors of the real 
output from the ideal output, and obtain a quick and accuracy 
respond for the power source that used in the particle acceler- 
ators. 

Unlike conventional control theories, MPC can forecast the 
output according to information of both the future reference 
input and the last output. By optimally determining the future 
input on every sampling period, the proposed method can ad- 
just the input in advance to diminish influence of the unex- 
pected factors as much as possible, so that the real output can 
trace the ideal output closely. The control strategy employs 
the state equation as the predictive model, and noises of dif- 
ferent frequency are added to the system to search the impact 
of the MPC on the system. The proposed control scheme is 
simple and computationally efficient, requiring just a few op- 
erations. The results of simulations and experiments show 
improvement on reduction of the noise influence. 
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